Phase transitions in interacting domain-wall model 



Jae Dong Noh and Doochul Kim 

Department of Physics and Center for Theoretical Physics, Seoul National University, 

Seoul 151-742, Korea 
(February 1, 2008) 

Abstract 

We investigate the interacting domain- wall model derived from the triangular- 
lattice antiferromagnetic Ising model with two next-nearest-neighbor inter- 
actions. The system has commensurate phases with a domain-wall density 
q = 2/3 as well as that of q = when the interaction is repulsive. The 
q = 2/3 commensurate phase is separated from the incommensurate phase 
through the Kosterlitz-Thouless (KT) transition. The critical interaction 
strength and the nature of the KT phase transition are studied by the Monte 
Carlo simulations and numerical transfer-matrix calculations. For strongly at- 
tractive interaction, the system undergoes a first-order phase transition from 
the q = commensurate phase to the incommensurate phase with q ^ 0. 
The incommensurate phase is a critical phase which is in the Gaussian model 
universality class. The effective Gaussian coupling constant is calculated as a 
function of interaction parameters from the finite-size scaling of the transfer 
matrix spectra . 
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I. INTRODUCTION 



The triangular-lattice antiferromagnetic Ising model (TAFIM) displays rich critical phe- 
nomena. The TAFIM with nearest-neighbor coupling i^T (< 0) is described by the Hamilto- 



nian 



-^ = '<Z.,s, (1) 

where (ij) denotes pairs of nearest-neighbor sites of triangular lattice and Sj = ±1. The 
ground states of the model are infinitely degenerate due to frustration on each elementary 
triangle and correspond to a critical state with algebraic decay of correlations Each 
ground-state configuration can be mapped to a state of the 'triangular solid-on-solid (TI- 
SOS)' model which describes equilibrium shape of a simple-cubic crystal near its [1,1,1] 
corner or growing of the simple-cubic crystal along the [1, 1, 1] direction 0], and also of the 
domain- wall model |1|,^. The latter describes the commensurate- incommensurate (C-IC) 
phase transitions [Q. When a two-dimensional system with anisotropic interactions has 
degenerate ground states, excitations from a ground state may take the form of domain 
walls of striped shapes. A typical example is the axial next-nearest neighbor Ising (ANNNI) 
model 0. Though domain- wall type excitations are energetically unfavorable, they also 
carry a finite amount of entropy. So the system is in an ordered phase without domain-wall 
excitation (C phase) at low temperatures while it is in a modulated phase with finite density 
of domain wall for sufficiently high temperatures. 

Ground-state properties of the TAFIM are studied by taking the limit K —>■ —oo in 
Eq. and by assigning statistical weights to each ground-state configuration. The new 
effective Hamiltonian of interest here can then be written as 

n = - J2 ^aSiSj - eaSiSj -hj^^i ■ (2) 
(iJ) ^ 

Here, (ij) {{{ij))) denotes pairs of nearest-neighbor (next-nearest-neighbor) sites, 6a (ca), 
a = 1,2,3, are direction-dependent nearest-neighbor (next-nearest-neighbor) couplings and 
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h is a magnetic field. Each allowed spin configuration {sj} is understood to be one of 
the ground-state configurations implicitly. The statistical weight for each configuration is 
proportional to e^''^^''^*'^'. 

If there are only nearest-neighbor couplings, the system is equivalent to the non- 
interacting domain- wall model and can be solved exactly |]T],^. The system exhibits 
Pokrovsky-Talapov (PT) transition which separates an ordered commensurate (C) phase 
and a disordered incommensurate (IC) phase. The IC phase is a critical phase which is in 
the Gaussian model universality class. The free-energy functional of the Gaussian model is 
given by 



(3) 



2 \dxi J 2 \dx2 J 

where Ki and K2 are the anisotropic stiffness constants. The Gaussian model has scaling 
dimensions 

^.,-| + f,^ (4) 



where p and q are integers and g = 2'K\JK\K2 is called the Gaussian coupling constant 0. 
The Gaussian coupling constant is equal to 1/2 throughout the IC phase in this limit 

Effect of the magnetic field for the system with = Ca = has been studied by numerical 
transfer-matrix calculations and Monte Carlo simulations 0]. The system undergoes 
Kosterlitz-Thouless (KT) phase transition from the disordered critical phase to the 



ordered phase as h increases. 

When there are two next-nearest-neighbor couplings (ei 7^ 0, £2 7^ 0, €3 = 0) as well as 
nearest-neighbor couplings, li, in Eq. (0) can be interpreted as a Hamiltonian for an inter- 
acting domain- wall system [^]. The nearest-neighbor and next-nearest-neighbor couplings 
play the role of the fugacity of domain walls and the interaction between domain walls, re- 
spectively. There are two types of domain walls and the two next-nearest-neighbor couplings 
ei and €2 control wall- wall interactions between respective types. In the previous paper [Q, 
we considered the partially interacting domain-wall (PIDW) model where domain walls of 



only one type interact with each other (ei = 0, €2 7^ case). The PIDW model is equivalent 
to the general five- vertex model which in turn is exactly solvable through the Bethe Ansatz 
method. We obtained exact phase diagram and quantitative understanding of the effect of 
the interaction in the critical properties of the PIDW model. When the repulsive interaction 
is strong, the system is in a novel C phase with domain- wall density g = 1/2 which is absent 
in the non- interacting system. In fact, this phase is equivalent to the C phase appearing in 
the ANNNI model where the interactions between domain walls are assumed to be infinitely 
repulsive The IC phase is also the critical phase in the Gaussian model universality class 
but the Gaussian coupling constant g varies continuously throughout the IC phase. 

In this paper, we study the interacting domain-wall (IDW) model where domain walls of 
both types interact with same strength (ei = 62 7^ case). We are concerned with the nature 
of the ordering and the phase transition. We find that the system has an ordered C phase 
with g = 2/3 when the repulsive interaction is sufficiently strong and suggest that the phase 
transition between the IC phase and the g = 2/3 C phase is the KT transition. We argue 
that the transition is described in the continuum limit by the Gaussian model Eq. (^) with 
a symmetry-breaking field V3 cos (30(r)). The repulsive interaction acts as the symmetry- 
breaking field. When V3 is sufficiently small, the system is in the critically disordered state 
with scaling dimensions given in Eq. @). The Gaussian coupling constant is not constant 
any more in the critical region but varies continuously since it is renormalized by the V3 field. 
The symmetry-breaking field V3 has a scaling dimension X3 0. As V3 increases, it becomes 
relevant in the renormalization group sense. This happens when X3 = 2 (or g = 9/4). 
The phase transition involved in this order-disorder transition is the KT transition. So, we 
predict that the KT transition occur as the repulsive-interaction strength varies. 

On the other side, where the attractive interaction is strong, the domain walls are 
bounded with each other. We find that, as the chemical potential of domain walls increases, 
there is a first-order phase transitions from the g = C phase to the IC phase accompanied 
by a finite jump in g. 

This paper is organized as follows. In Sec. || we introduce the IDW model and present 
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the relation between it and the TISOS modeL Possible ordered phases in the IDW model 
are discussed. Unlike the PIDW model, the IDW model is not exactly solvable. So, we 
investigate the IDW model by numerical diagonalization of the transfer matrix and Monte 
Carlo simulations. In Sec. ^ we analyze the eigenvalue spectra of the transfer matrix 
using the finite-size scaling theory. The transfer matrix applies to a system with a cylinder 
geometry: infinitely long in one direction, and periodic with a finite size in the other 
direction. The calculations are performed for the strip width up to = 18, corresponding 
to a transfer matrix of size 2^^ x 2^^. Using the finite-size scaling of the eigenvalue spectra, 
we calculate the Gaussian coupling constant g from which the KT transition point for 
the repulsive interaction is located. The first-order phase transition line for the attractive 



interaction is obtained from the eigenvalue spectra and an analytic approximation. In Sec. |V 
results of Monte Carlo simulations are presented. The simulations are done for a lattice of 
size up to 120 X 120 while the domain- wall density is fixed to 2/3. As well as the specific 
heat, we also calculate the fluctuations of the position and the density of domain walls since 
they serve as order parameters for the KT transition. Monte Carlo results are consistent 
with those of transfer-matrix calculations. In Sec. |V| we summarize our results and discuss 
the origin of the difference between the PIDW model and the IDW model. 



II. THE MODEL AND ITS GROUND STATES 

We consider two-dimensional domain-wall model on a finite N x M triangular lattice 
with one side of elementary triangle lying along the horizontal direction. The configurations 
of domain walls are mapped from the ground-state configurations of the TAFIM [^]. The 
domain walls propagate vertically along the two directions of elementary triangle. We call 
the domain wall of type 1 (2) if it is left (right) moving. A single chain of domain wall 
has an entropy log 2 per unit length since there are two ways for domain walls to move. 
The configurations of domain walls are subject to the restriction that they never be created 
or annihilated during propagation. Thus, we do not allow the 'dislocations' which form 
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topological excitations. However, effect of such excitations will be discussed briefly in Sec. 0. 

Domain walls have the chemical potential fi and interaction energy e per unit length. 
Relationships between {fi,e) and (5a, Ca) of Eq. are derived in Ref. ||^. The interaction 
energy is assigned to each segment of adjacent-parallel domain walls. The grand-canonical 
partition function of the IDW model is written as 

N 

^G(x,y) = E^*'^ E (5) 

Q=0 {pi,...,pq} 

where Q is the number of domain walls per row, x = e'^ is the fugacity of a domain wall 
segment of unit length and y = is the interaction parameter. Here, we use the convention 
that the 'temperature' of the domain-wall system is absorbed into /i and e. Note that this 
temperature is different from that of the TAFIM. The function L(pi, . . . counts total 
length of adjacent-parallel domain walls where Pi{j) denotes the horizontal position of the 
i-th domain wall at vertical position j. Since configurations with different number of domain 
walls do not couple with each other at all, the partition sum can be evaluated independently 
for each Q] 

N 
Q=0 

where Zc is the canonical partition function, representing the second sum in Eq. (^), and 
q = Q/N is the do main- wall density. The grand-canonical and canonical ensembles are 
connected through the relation 

1 d log Zc , . 

^ = -MN^^ ■ 

Before proceeding the analysis of the IDW model, we digress to explain the equivalence 

of the domain-wall model and the TISOS model. The TISOS model is a solid-on-solid model 

on a triangular lattice on which integer-valued heights hi at site i = (xj, yi) are assigned 

with the restriction that the height differences between any nearest-neighbor sites are ±1 or 

±2. See Ref. ll| for details. Let be the number of domain walls on the left-hand side of 

the point {xi,yi). Then, the mapping 

hi = 3ni-2xi (7) 



relates a domain-wall configuration to a height configuration. Note that the height variables 
satisfy the periodic boundary condition /^(xi+Ar.j/i) = only if the domain-wall density 

g = 2/3. In this case the TISOS model so constructed, describes the equilibrium shape of 
a cubic crystal viewed from the [1, 1, 1] direction. The case of g 7^ 2/3 corresponds to that 
viewed from other directions. 

'Ground states' of the model for given x and y are the configurations which maximize the 
Boltzmann weight in Eq. (|^). For the attractive interaction (e < 0), a C phase with q = 
(1) is the ground state when fi < e > e). There will be a first-order phase transition 
along the line /i = e at zero temperature, i.e., \e\ — > 00. On the other hand, for the repulsive 
interaction (e > 0), q = 0(1) C phase is the ground state when /i < 0(/U > 3£:). When 
< fi < 3e, the g = 2/3 C configurations become the ground states. These new ground 
states correspond to the most densely-packed domain-wall states without interacting pairs 
and are three-fold degenerate. They are shown in Fig. 1. 

To understand finite-temperature behaviors, we should consider excitations from these 
ground states. Excited states from the g = state are formed by domain-wall formations. 
If the domain walls are far apart from each other (this is the case when the attractive 
interaction is not so strong), the energy cost is — /i and the wandering entropy gain is log 2 
per unit length of domain walls. Thus there will be a phase transition at /i = — log 2. 

An excited state of the g = 1 state is formed by removing one segment of domain walls 
from each row. We will call it the domain-wall hole excitation. Unlike a domain wall, the 
domain-wall hole need not form a continuous path. Thus a single domain-wall hole has a 
wandering entropy logA^ and takes energy cost fi per unit length. So the g = 1 state is 
always unstable under domain-wall hole excitations at any value of fi. 

The g = 2/3 C state has complicated excitations. The simplest kind of excitations is 
generated by locally reversing two units of a domain wall with respect to a vertical axis. 
Its excitation energy is 26. This excitation is a local excitation so it cannot destroy the 
long-range order at low temperatures. Another important excitation is formed by domain 
boundaries between the three ground-state configurations. This is illustrated in Fig. 2 where 
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three ground-state configurations of Fig. 1 coexist forming domain boundaries. 

Nature of such excitations becomes more transparent when we introduce a coarse-grained 
field variable as follows. First we assign a height variable hi to each lattice site according to 
Eq. d^). The height variables corresponding to the domain wall configuration of Fig. 2 are 
also shown in Fig. 2 by small numbers. We now group lattice sites into a triangular array of 
elementary triangles shown as dotted triangles in Fig. 2. To each elementary dotted triangle 
a, we then assign a coarse-grained height variable 



ha — {ha\ ~\~ ha2 ~\~ ho,^)/3 



(8) 



where a^'s denote the three sites of the triangle a. The values of h^ for the configuration 
of Fig. 2 are shown by big numbers in the middle of each dotted triangle, ha take integer 
values and ha = ha' (mod 3) if local configurations at a and a' are in the same ground-state 
configuration. At each boundary, Aha = ±1 or ± 2. In the continuum limit, the critical 
behaviors of the system are described by the Gaussian model with free-energy functional JF 
in Eq. (^ if there is no interaction between domain walls. The repulsive interaction favors 
the coarse-grained height h{r) to be an integer where h{r) denotes ha in the continuum 
limit. So we postulate that the order- disorder phase transition in the IDW model with 
q = 2/3 is described by the Gaussian model with a symmetry-breaking field whose free- 
energy functional is given by 



2^^ [dxi 



+ V3COS (3(/)(r)) 



(9) 



where 0(r) = -^h(r) and the parameters Ki,K2 and V3 are functions of the interaction 
strength y. It is known that this system has two phases: critically disordered phase (high 
temperature phase) and the ordered phase (low temperature phase) |T0|. In the critical 
region the model is renormalized to 



2 \dxi 



(10) 



where Ki and K2 are renormalized stiffness constants which are functions of Ki, K2 and V3. 
The scaling dimension of the symmetry-breaking field is X^^ = 9/{2g) with the effective 
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Gaussian coupling constant g = 2t:\I K1K2. In the critical region the symmetry-breaking 
field is irrelevant (i.e., Xs^o > 2). The phase transition into the ordered phase occurs at 
the values of the model parameters which satisfy the relation X3 = 2 and belongs to the 
universality class of the KT transition. The transition point can be found from the condition 
9 = 9/4. 

In this section, we presented ground states of the IDW model and the conjecture that 
the C-IC transition of the IDW model with the repulsive interaction at g = 2/3 is described 
by the free-energy functional given in Eq. (^. In the next two sections, we will confirm our 
conjecture by transfer matrix calculations and Monte Carlo simulations. 

III. FINITE-SIZE SCALING ANALYSIS OF TRANSFER-MATRIX SPECTRA 

We study the IDW model on a finite lattice with width N and height M using the 
transfer matrix. The grand-canonical partition function Zq in Eq. (j^) is written as 

N 

Zg=Y1 ^"""^ Tr Tj^ 

Q=0 

where Tg is a transfer matrix which acts on the nCq dimensional space spanned by statekets 
|ni, . . . , hq) specifying positions of Q domain walls. The matrix element (n'^, . . . ^u'q \ Tg | 
77.1, .. . ^uq) represents the Boltzmann weight for a domain- wall configuration between two 
rows of a lattice as shown in Fig. 3. 

We denote the p-th largest eigenvalue of Tg by e~^Q and Eq will be called the energy 
of the ^th excited level. In the limit M — > cxd, the canonical partition function is given by 
the largest eigenvalue, i.e., Zc = e"*^'^^. In the IC phase, q = Q/N varies continuously as 
a function of /i and e. When q is given, the chemical potential should be chosen so as to 
satisfy the finite-size version of Eq. (H) 

^~ 2N 

The eigenvalue spectrum for a conformally invariant system follows the universal scaling 
form 
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Re E^N) = NU + ^ (X^ _ o f 1 ) (11) 



N \ 12/ VA^. 
where c is the central charge, is the scahng dimension of the operator associated with 



eigenstate a, ( is an anisotropy factor, and f^o is the bulk free energy per site |TT[. The IC 
phase of the domain-wall system is the critical phase in the universality class of the Gaussian 
model whose central charge is 1 and scaling dimensions are given in Eq. In the previous 
paper we showed that the p-th excited level of Tq leads to the scaling dimension Xp o and 
excitations of the largest eigenvalues of Tq^g relative to that of Tq lead to Xq ^. Using this 
knowledge, we can estimate the Gaussian coupling constant by an estimator g{N) defined 
by 

9iN) 



Re(E0+i + E0_i-2E0; 



\ 2Re(E^-E0) 
which should approach g as N oo. A sequence of g{N) is obtained by numerical diago- 

nalization of transfer matrix for a system of width N up to 18. We present the estimates in 
Fig. 4(a) and Fig. 4(b) for the repulsive {y < 1) and attractive {y > 1) case, respectively for 
several values of y. The curve for y = 0.2 in Fig .4(a) is noticeable. g{N) apparently does 
not converge to a finite value as N increases when q = 2/3. It means that when y = 0.2 
the eigenvalue spectrum does not follow the finite-size scaling form of Eq. (0) and that the 
system is no longer in a critical phase. We gave an argument in the previous section that the 
system undergoes a phase transition into the ordered phase as the repulsive interaction be- 
comes strong and that the Gaussian coupling constant should be 9/4 at the transition point. 
Thus we estimate the critical interaction strength y^ by solving the equation g{N) = 9/4 
numerically. For this purpose, we find that the alternative estimators defined by 

^ 6(g^(iV + 3)-E°(iV)) 

(i5;^(iV + 3)-EO(iV + 3))-(E^(iV)-EO(iV)) 

converge better. The solutions of g{N) = 9/4, denoted by yi{N), are given in Table I. The 
strip-width dependence of yi{N) follows from corrections to scaling of the eigenvalues. At 
the KT transition point there are logarithmic corrections due to the presence of a marginally 
irrelevant field Thus the sequence of estimator yi{N) is extrapolated by fitting it to the 
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form yi{N) — 1/2 + a/{b + logN). Since the sequence 1/2 thus obtained also show slow- 
convergence, we fit it again to the same form and obtain the result yc — 0.252 ± (0.003). See 
Table I. 

When y < y^ there appears the q — 2/3 C phase. As in the q — 1/2 C phase of the 
PIDW model, q remains locked to the value 2/3 in the range /x_ < < Estimates of 
the chemical potential at the lower and upper boundaries of the C phase are given by 

where Qo = 2N/3. Since they are evaluated from finite-size systems, we should know the 
finite-size scaling form of iJ,±{N). Due to the marginal corrections at the phase boundary, 
the largest eigenvalues ior Q — Qo ± m is expected to have scahng forms 

It then follows that iJi±{N) behaves as 

Accordingly, we estimate /x± by fitting ij,±{N) to the equation IJ,±{N) = iJ,±+a/N+b/N log N. 
The result is shown in Fig 5. 

We have thus determined the phase boundary of the q = 2/3 phase using the finite-size 
scaling of the eigenvalue spectra of the transfer matrix. Next, we consider the system near 
the g = C phase. As the chemical potential increases, domain walls begin to form. They 
form a free or bound state depending on the interaction strength y. To see this, we calculate 
in the Appendix the ground-state energy and un-normalized eigenket of Tq for Q = 2 from 
the Bethe Ansatz method. The results for the maximum eigenvalue and the eigenket are 

= -21og2, \X)^j:m<n\m,n) V < I ^^^^ 

El = - log {2y + ^) , I A) = E^<n (2(y - l))-^"-'") |m, n) if y > | , 
where the thermodynamic limit is taken. The eigenket has a qualitatively different properties 
on both side of y = 3/2 hue. 
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As can be seen from Eq. ([T2D, domain walls form a free state when y < 3/2. So, the 
C-IC transition is the PT transition as in the non-interacting and the PIDW models. The 
domain wall density has a square-root dependence on the chemical potential 



g ~ (/i - /i 



1/2 



with the critical chemical potential /ic = — log 2. 

When y > 3/2, Eq. (|T2|) shows that domain walls form a bounded state. The two domain 
walls are bounded with the mean distance 

The ground-state energy of the bounded state consists of two parts E2 = -Ec.m. + -E'int. where 
Ec.m. = — log 2 comes from free motions of the center of mass of the two domain walls and 
Eint. = log {y + l/4(y — 1)) from the effective interaction. To generalize Eq. ([T^ ) for Q > 3 
where exact results are not available, we construct a simple approximation for Eq as a sum 
of Ec.m. and -Eint.- Thus we write the ground-state energy of the Tq for Q 7^ 2 as 

E'q = ^c.m. + (Q - l)Emt. (13) 

1 \ 



log2-(g-l)logh/ 



4(1/ -1)7 

This approximation will be called the two-domain-wall approximation. It should be valid 
as long as q « l/do. When q > 1/do, domain walls feel additional statistical repulsive 
interactions. This statistical repulsions then drive the system into the IC phase. Within 
the two-domain-wall approximation, the domain-wall system undergoes a first order C-IC 
transition from the g = C phase to the IC phase. Since Eq. (0) is linear in Q, the 
transition occurs at the chemical potential 

''« = -i°i5(y + j^) (14) 

while go; the discontinuity of q at the transition, can be taken to be 1/do approximately. 

/io and go are also determined numerically. When y > 3/2, Eq for finite obtained from 
numerical diagonalization of Tq becomes concave downward in some region < Q < Qq 
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indicating a first order phase transition. We took the convex envelop of Eq and estimated 
yUo and go from the slope of the straight line enveloping the concave region and from the 
point where the straight line and the curve of Eq are tangential, respectively. The results 
are shown in Figs. 6 for five values of A^. These are compared with Eq. ( [T^ ) and the 
approximation 1/do in Figs. 6. 

Resulting phase diagram of the IDW model as obtained from the transfer-matrix calcu- 
lations is summarized in Fig. 7. This is the main result of this work. 



IV. RESULTS OF MONTE CARLO SIMULATIONS 

In the previous section, the assumption that the repulsive interaction induces the KT 
transition at the g = 2/3 C phase boundary is used to estimate the critical interaction 
strength t/c from the numerical diagonalization of the transfer matrix. In this section, we 
will confirm the nature of the phase transition more directly by Monte Carlo simulations at 
g = 2/3 on a finite L x L lattice. 

At the KT transition point, there is no diverging peak in the specific heat. So, we should 
study another quantities such as the renormalized stiffness constants Ki and K2 in Eq. (|10|). 
They can be evaluated from the linear response 

\dxk/ Kk 
produced by an additional free energy term 




where A; = 1 or 2. From the fluctuation-dissipation theorem, we see that 

90(r)90(r')\ /90(r)\ /90(r')\ 
dxk dxk / g \ dxk /q \ dxk / 
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where V is the volume and (■ ■ ■)o denotes the average with Vk = 0. These quantities are called 
helicity moduli and serve as order parameters for systems exhibiting the KT transition, e.g., 
XY model [0. 

Using Eqs. (0) and which connect the Gaussian model and the domain-wall model, 
Kk can be written in terms of the quantities of the domain-wall model as 

l_47r2 r//. 2.\/. 2.\\ /. 2.\/. 2 



Afc^j - [Aknj - -4,1 jy^ - y^fc^i - 3*^^'!/^ {^^knj - -4,i 



where A^nj = rii^ek — "^i, is the unit vector along the k direction and 4,1 is the Kronecker 
delta function. From the definition of rij, we see that 



AiUi 



1 , if there is a domain wall between (xj, i/i) and (xj, i/i) + ei 
, otherwise 



1 , if there is a domain wall of type 1 between (x,, i/i) and (xj, yi) + 62 
A2^i = { —I , if there is a domain wall of type 2 between (x,, yi) and (xj, yi) + 62 
, otherwise 

Using this relation, we can write the expression for the stiffness constants as 
1 47r2 



1 47r2 //Q2-gi\'\ 47r2 



EAX. (15) 



where Q {Qi,2) is the total length of the domain walls (of type 1 or 2) and AXi = pi{L)—pi{0) 
is the difference of the position of the i-th domain wall at the top to that at the bottom. 
The equations show that Ki controls the fluctuations of the domain-wall length and K2 the 
fluctuations of the domain-wall positions. 

We performed simulations with the periodic boundary condition along the horizontal 
direction and the free boundary condition along the vertical direction with domain-wall 
density fixed to 2/3. Updating rule from a domain- wall configuration {pi{j)} to {p'i{j)} is 
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as follows. There are four kinds of configurations for any 

it) piij - 1) + 1/2 = piij) = piij + 1) - 1/2 

(ii) pi{3 - 1) - 1/2 = p,{j) = pi{j + 1) + 1/2 

(m) - 1) - 1/2 = = + 1) - 1/2 

{iv) pi{j - 1) + 1/2 = pi{j) = pi{j + 1) + 1/2 . 

A given configuration can be changed locally only for the cases (in) and (iv). For the 
case (m), Pi{j) is updated to = Pi{j) + 1 with a probability provided 

7^ — 1, while for the case {iv), Pi{j) is updated to p[{j) — Pi{j) — 1 with a 
probability peovided Pi-i{j) Pi{j) — 1- For other cases, the configuration is 

not altered. Noting that the updating trial at has nothing to do with domain walls 

(i ± 2, J ± 2), we divide the set into four parts depending on whether i,j are even 

or odd and choose randomly one of the four subsets and update the whole domain walls in 
that subset. 

The free boundary condition along the vertical direction is necessary to calculate the 
quantities AXj. On the other hand, one should adopt, in principle, a grand canonical 
ensemble in which the domain-wall number could vary in other to calculate the fluctuations of 
the domain-wall number. But it poses too much technical difficulty, so we measure instead 
fluctuations of the domain- wall number in a fixed rectangular region of size {L/2) x L 
during simulations of the L x L system. For such an ensemble, an elementary probability 
consideration shows that 1/Ki should given by 

^ = ^ m - («r) 

where Q' is the number of domain walls inside the {L/2) x L region. 

In Fig. 8, we present results for the energy fiuctuations C which is defined by 

c ^ ^ imp}) - iuip}))?) 

for systems of linear size L = 12, 36, 60, 90 and 120. As expected, there is a converging peak 
instead of a diverging peak. The peak position is below the estimated critical interaction 
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strength, i.e., at the low-temperature region. This is consistent with the behaviors in the 
XY model where the peak is at the high-temperature region. The KT transition in the 
XY model is the vortex-antivortex unbinding transition and that in the IDW model is a 
spontaneous symmetry-breaking transition. The effects of the vortex-antivortex and the 



symmetry-breaking field have a dual relation WU 



where K is the inverse temperature, yo is the fugacity of the vortex and antivortex, and yp 
is the strength of the symmetry-breaking field. We see that the high and low temperature 
regions are inverted. 

Figs. 9(a) and (b) show the measured values of l/i^i and I/K2, respectively, as a function 
of y for several values of L. Though we did not try to extract the quantitative results, 
they show manifest crossover behaviors near y ~ yc- They become smaller and smaller in 
the ordered regions as the system size L increases and the points at which the crossover 
behaviors set in become closer and closer to the estimated critical-interaction strength y^. 
These behaviors are in accord with Monte Carlo results of the XY model |T2|. Finally, 



Fig. 9(c) shows the values of g obtained from the relation g = 271^X1X2. The coupling 
constant estimated in this way is in good agreement with that obtained from the transfer 
matrix spectra. That the values of g at y^. fall short of the expected value 9/4 is attributed to 
slow convergence arising from the presence of logarithmic corrections. From these numerical 
results, we conclude that the phase transition belongs to the universality class of the KT 
transition. 



V. SUMMARY AND DISCUSSION 

In summary, we have investigated the phase transitions in the IDW model. There are 
three phases: the g = and g = 2/3 C phases and the IC phase. The transition between 
the g = C phase and the IC phase is the PT transition when the interactions between 
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domain walls are repulsive or weakly attractive. When the interaction is strongly attractive, 
it becomes the first-order phase transition with a discontinuity of the domain-wall density. 
The nature of the transitions are altered because of formation of bounded domain-wall 
states. As to the transition between the g = 2/3 C phase and the IC phase, we suggested 
the free-energy functional is that of the Gaussian model with a symmetry-breaking field 
with the scaling dimension X^^q. The symmetry-breaking field accounts for effects of the 
repulsive interaction and induces the KT transition from the critical phase to the ordered 
phase. We confirmed this scenario by calculating the stiffness constants from the Monte 
Carlo simulations. From the transfer matrix method, we calculated the Gaussian coupling 
constant in the IC phase and the critical interaction strength. Phase boundaries in the e-n 
plane are also determined from the transfer matrix method. 

The domain-wall model is constructed from the TAFIM with the restriction that three 
spins on elementary triangles cannot have the same sign. If this is allowed the excitation 
of this type constitutes a topological excitation where two-domain walls are created or 
annihilated [Q. It plays the role of vortex and antivortex excitation in the XY model. 
The scaling dimension of the excitation is Xq^2- Since it is irrelevant when g > 1, it cannot 
alter the nature of the C-IC transition of the g = 2/3 C phase where g = 9/A pTSf . 

In the PIDW model where the interaction and the chemical potential is anisotropic 
the repulsive interaction stabilizes the q = 1/2 C phase. But, the phase transition is of 
the FT type. Thus the anisotropy plays an important role in determining the nature of the 
phase transition. The anisotropy produces an extra term in the free-energy functional which 
couples to the differences of the number of domain walls of type 1 and 2. So the FIDW 
model would be described by the effective free-energy functional 



where the symmetry-breaking field accounts for the ordering induced from domain-wall 
interactions and the linear term —pd(j)/dx2 accounts for the effect of the anisotropy. In this 
model, the phase transition is induced by the coupling p. It induces a domain-wall type 
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excitation and the phase transition is in the universahty class of the PT transition [Q. 

Even though we studied only the equilibrium properties of the IDW model in this work, 
non-equilibrium version is also of interest since it describes a driven domain-wall model or 
the hypercube stacking model 0] and also the charge conduction in the charge-density-wave 
system []T^. We leave it for further study. 
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APPENDIX A: 

In this appendix, we derive expressions for the largest eigenvalue A and the corresponding 
eigenket |A) of the transfer matrix Tq for Q = 2 analytically. The eigenket is of the form 

|A) = ^ a{m,n)\m^n) 

l<m<n<N 

where \m, n) denotes a state with two domain walls at sites m and n. The eigenvalue 
equation T2IA) = A|A) gives linear equations for coefficients a{m,n): 

Xa{m,n) = ^ a{m' ,n'){m,n\T2\m' ,n') . (Al) 

m'<n' 

When T2 is applied to a stateket \m,n), it generates a new stateket: 

|m, n) + |m, n + 1) + |m + 1, n) + |m + 1, n + 1) if n > m + 1 
y\m, n) + \m, n + 1) -\- y\m + 1 , n + 1) if n = m + 1 

where y is the interaction parameter. Thus Eq. becomes 



T2|m, n) 



Aa(m, n ^ m + 1) = a{m, n) + a{m — 1, n — 1) + a{m, n — 1) + a{m — 1, n) (A2) 
Aa(m, n = m + 1) = ya{m, m + 1) + ya{m — 1, m) + a(m — 1, m + 1) . (A3) 

These equations are solved by the Bethe Ansatz 

a(m, n) = AuzT""z2'' + ^2i^r"^2"'" (A4) 

where y4i2, A21, Z\ and ^2 are constants to be found. Eq. (|A2|) is automatically satisfied by 
a(m, n) given in Eq. ( |A^ ) provided the eigenvalue is 

A= (1 + Zi)(l + Z2) . 

We can make Eq. ([A3| ) be satisfied by imposing the condition that 



[a(m, n) + a(m — 1, — 1) + a(m, n — 1) + a(m — 1, 

= ya{m, m + 1) + ya{m — 1, m) + a(m — 1, m + 1) 

19 



This gives the relation between A12 and A21 as 

Au l-(y-l)(ir + ^2 



A21 + 

If we use a periodic boundary condition a(n, m + N) = a{m, n), we obtain the equations for 



zi and Z2: 



N _ A21 TV _ ^12 

^12 ^21 



Multiplying the two equations together gives {ziZ2)'^ = 1 which implies that Z1Z2 = t where 
T is an A^-th root of unity. The eigenvector corresponding to the maximum eigenvalue comes 
from the solution with r = 1 which means that zi = 1/ Z2 = z. The value of z is determined 
from the equation 

z 

The solutions of Eq. ([A5|) is given by 



{2y-2)^ ) " y ^ 2 

i + ify<| . 

In the limit — > 00, this gives the eigenvalues and eigenvectors presented in Eq. (|T^). 

Unfortunately the Bethe Ansatz works only for Q = 2 and fails when Q > 3. Thus we 
are not able to get analytic results for general Q. However, the solution for Q = 2 provides 
a useful approximation near the g = C to IC transition as discussed in Sec. |ITT|. 
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Table Caption 



Table I. The transfer-matrix calculation of the KT transition point of the IDW model. 

Estimates of the critical interaction strength yi{N) is found by solving the equation g{N) = 
9/4 by the Newton iteration method. The sequence of yi{N) is fitted to the form yi{N) — 
y2 + a/{b + logN) to obtain a sequence of extrapolated values y2{N). yc is estimated by 
fitting y2{N) again to the same asymptotic form. The error is estimated from the difference 
between y^ and 1/2(9) 
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Figure Captions 

Fig. 1. Three ground states of the domain-wall system with repulsive interactions when 
< n < 3e. The broken line represents a triangular lattice and solid line represents sections 
of the domain wall. Each ground state is obtained by periodic repetitions of each unit in 
the vertical direction. 

Fig. 2. An excited state of the domain-wall system with repulsive interactions. The small 
numbers on vertices of a triangular lattice denote heights hi of an equivalent TISOS config- 
uration and the large numbers in dotted triangles denote the coarse-grained heights ha- 
Fig. 3. An example of a domain- wall configuration in successive two rows of a triangular 
lattice of width A'^ = 6. The domain- wall configuration of the lower row is represented 
by the stateket |1,2,3,5) while that of the upper one is represented by |1,3, 4,5). The 
transfer-matrix element (1, 3, 4, 5|Tq=4|1, 2, 3, 5) denotes the Boltzmann weight of the above 
configuration and is equal to y. 

Fig. 4(a). Estimates of the Gaussian coupling constant g{N) as a function of q for several 
values of y. The symbols A, O, +, □, and • correspond to A^ = 6, 9, 12, 15, and 18, 
respectively. The values of y are 0.2, 0.4, 0.6, 0.8, and 1.0 from top to bottom at each values 
of q. The broken fine denotes the line g — 9/4:. (b) Same as in (a) but for the values of y 
1.0, 1.1, 1.2, 1.3, 1.4, and 1.5 from top to bottom. 

Fig. 5. Phase boundary of the q — 2/3 C-IC transition in the /i-y plane. /i± are obtained by 
extrapolating ijl±{N). The broken line is drawn at the critical interaction strength y — 0.252. 
fi+ and fj,- merge into a single value near that line. Solid lines are guides to the eye. 
Fig. 6. Estimates of Hq (a) and qo (b) obtained from transfer- matrix calculations (symbols). 
The two-domain- wall approximation for /iq and the approximation qg l/dg are also shown 
in (a) and (b), respectively, by solid lines. 

Fig. 7. Phase diagrams of the IDW model with the repulsive {s > 0) and attractive {s < 0) 

interactions. The solid lines ( — ) denote the PT transition, the broken line ( ) the KT 

transition and the dotted fine (■ ■ ■) the first-order transition. The PT transition lines are 
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straight lines with slope — l/log2. 

Fig. 8. Monte Carlo results for the energy fluctuation C. This result indicates that the 
energy fluctuation remains flnite when the system enters the ordered phase. The broken line 
is drawn at y = 0.252. 

Fig. 9. Monte Carlo results for the renormalized stiffness constants Ki (a), K2 (b) and the 
Gaussian coupling constant g (c). They show characteristic behaviors of the KT transition. 
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TABLES 



N yi 2/2 Vc 

3 0.170790486223 0.265782 0.252 

6 0.197135051080 0.256660 

9 0.206717177287 0.255360 

12 0.211836668182 

15 0.215121556318 



Table I 
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Fig. 9(a) 




Fig. 9(b) 




Fig. 9(c) 



